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ABSTRACT 



Context. Whether the Cygnus X complex consists of one physically connected region of star formation or of multiple independent 
regions projected close together on the sky has been debated for decades. The main reason for this puzzling scenario is the lack of 
trustworthy distance measurements. 

Aims. We aim to understand the structure and dynamics of the star-forming regions toward Cygnus X by accurate distance and proper 
motion measurements. 

Methods. To measure trigonometric parallaxes, we observed 6.7 GHz methanol and 22 GHz water masers with the European VLBI 
Network and the Very Long Baseline Array. 

Results. We measured the trigonometric parallaxes and proper motions of five massive star-forming regions toward the Cygnus X 
complex and report the following distances within a 10% accuracy: 1 .30*°-°'!, kpc for W 75N, 1 .46+°* kpc for DR 20, 1 .50+5J j5^ kpc 
for DR21, 1.36!j; |jkpc for IRAS 20290-1-4052, and 3.33!j;H kpc for AFGL2591. While the distances of W75N, DR20, DR21, and 
IRAS 20290+4052 are consistent with a single distance of 1.40 ± 0.08 kpc for the Cygnus X complex, AFGL2591 is located at a 
much greater distance than previously assumed. The space velocities of the four star-forming regions in the Cygnus X complex do 
not suggest an expanding Stromgren sphere. 



Key words. Masers - Astrometry - Techniques: high angular resolution - Stars: 
individual objects: W75N, DR20, DR21, IRAS 20290+4052, and AFGL2591 



formation - ISM: kinematics and dynamics. 



, 1. Introduction 



In the early days of radio astronomy, a conspicuously strong, ex- 
tended source of radio emission was foun d at Galactic longitude 
~80° and named the Cygnus X region (iPiddington & MinnettI 
11952), which also stands out in infrared surveys of the Galaxy 
(Odenwald & Schwartz 1993; see e.g., the spectacular Spitzer 
imaging in Kumar et al. 2007). All phases of star formation and 
stellar evolution are observed projected across the Cygnus X 
region, including a population of dense, massive, and dusty 
cores with embedded protoclusters and high-mass protostellar 
objec ts (Sridharan et al. 2002; Beuther et al. 2002; Mott e et al 
2007|). ultracompac t HII regions (jP ownes & Rinehartj Il966t 



type stars (of which ~65 O-type; IWright et al.ll201(i| and refer- 
ences therein), and some supernova remnants ( Uvaniker et al.l 
I2001h . 

The proximity of a large number of OB associations and 
molecular cloud complexes on the sky motivated the explana- 
tion of Cygnus X as consisting of various objects at different 
distances seen superposed on one another (e. g., Dickelet _d] 
I1969L IPinenbrink & Wendkerill988l and lUvaniker et aLlbOoTT 



Wendker et aOl99lHCvganowski et al...2003.) . hundreds of OB- 



Recently, the CO imaging survey of ISchneider et al. ( 2006h has 
rekindled the idea of Cygnus X as one large star-formin g com- 
plex, which was akeady suggested in the sixties by, e.g. IVeronI 
(1965). Obviously, both scenarios depend strongly on the dis- 
tances measured to the individual parts of the Cygnus X com- 
plex, because sources at a galactic longitude of ~ 80° could be in 
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the Local Arm (also named the Orion or Local Spur) and nearby 
(~l-2 kpc), in the Perseus Arm at ~5 kpc, or even in the Outer 
Arm at distances of ~10 kpc (e.g., the Cygnus X-3 microquasar). 

Unfortunately, distances to the Cygnus X objects are very 
difficult to obtain and have large uncertainties. First, the 
Cygnus X OB associations are too far for a parallax mea- 
surement with the Hipparcos satellite, which measured dis- 
tances of nearby OB associations out to a distance of 650 pc 
dde Zeeuw et al.l IT999I) . Second, at the Galactic longitude of 
Cygnus X, the radial velocity difference between the Sun and 
the Cygnus X region is close to the typical velocity disper- 
sion of interstellar gas in a high-mass star-forming region (SFR) 
(l-7km s"', Moscadelli et al. 2002) for distances up to 4kpc. 
Therefore, kinematic distances, which depend on the radial ve- 
locity, are not reliable for distances below 4 kpc toward this 
longitude, and most distance measurements rely on the spec- 
troscopy and photometry of stars. But, these estimates are also 
affected by large uncertainties (> 30%) becaus e the extinction 
towar d Cygnus X is very high and variable (ISchneider et al.l 
I2OO6 I). While one can find distance estimates between 1.2 and 
2 kpc in the literature, the generally adopted value is 1.7 kpc, 
based on the spectroscopy and photometry of the sta rs in 
the Cyg OB2 association b y iTorres-Dodgen et all (Il99l i) and 
iMassev & Thompson! (|1991|). More recently, a nearer distance 
of 1 .5 kpc was obtained by " Ikansonl (l2003h using new MK stel- 
lar classification spectra of the stars in Cyg OB 2. Of course, 
using a Cyg OB 2 distance for the entire complex assumes that 
Cygnus X is a physically connected SFR. 

The distance to this "mini-starburst" region is crucial for the 
many star formation studies performed toward it for its rich- 
ness and short distance from the Sun. As said, the distances are 
very uncertain: for example, three OB associations (Cyg OB 1, 
8, and 9) have distance estimates between 1.2 and 1.7 kpc, a 
difference of more than 30%. This distance range of 500 pc is 
almost ten times more than the extent of the Cygnus X region 
on the sky - 4° by 2° or 100 by 50 pc at a distance of 1.5 kpc. 
Therefore, important physical parameters of objects in this re- 
gion, such as luminosities and masses are uncertain by a factor of 
~1.7 given a distance uncertainty of 30%. For some of the SFRs, 
the distance estimates have a much wider range, namely from 
~1 kpc (AFGL2591) to 3kpc (DR21). To distinguish whether 
all clouds are at the same distance or are only projected close 
together on the plane of the sky, a direct estimate of distances to 
distinct objects toward Cygnus X is required. 

In this context, we used strong 6.7 GHz methanol masers 
and 22 GHz water masers as astrometric targets to measure dis- 
tances of five distinct SFRs toward the Cygnus X complex. No 
previous parallax measurements were carried out toward this 
region. Using the terminology of I Schneider et al.l (l2006h . the 
Cygnus X complex is divided about the Cyg OB 2 cluster at 
{I, b) = (80.°22, -i-0.°80) into a northern region, at Galactic lon- 
gitudes greater than about 80°, and a southern region, at lower 
longitudes. Methanol maser emission was observed toward four 
SFRs with the European VLB! Network (EVN): W 75N, DR 20, 
and DR21 in Cygnus X North, and IRAS 20290+4052_( which 
is likely part of the Cyg O B 2 association, see lOdenwaldll 19891 
[Parthasarathy et all 1 19921 hereafter IRAS 20 290) in Cy g nus X 
South. The pioneering work presented in iRvgl et all (1201 Ol) 
demonstrates the capability of the EVN to achieve parallax ac- 
curacies as good as 22 yum. Water maser emission was observed 
toward AFGL2591, which is projected within Cygnus X South, 
with the Very Long Baseline Array (VLBA). 

The observations presented here comprise one year of VLBA 
observations and two years of EVN observations with the addi- 



tion of Japanese antennas (the Yamaguchi 32-m antenna and the 
Mizusawa station of the VLBI Exploration of Radio Astrometry, 
VERA) for several epochs to increase the angular resolution. 
We report on the preliminary distances from the first two-year 
results, using only the EVN antennas, and the evidence that 
AFGL2591 is really projected against the Cygnus X region, 
hence not part of it. To optimize the results with the long base- 
lines afforded by the Japanese antennas, we need to develop 
more sophisticated calibration procedures, since several maser 
spots resolve on the ~9000km baselines. This analysis will be 
included in a future paper once the observations are completed. 



2. Observations and data reduction 

2.1. EVN methanol maser observations 

The EVrvQ observations were carried out in eight epochs be- 
tween March 2009 and November 2010 under project EB039. 
These dates were scheduled near the minima and maxima of the 
sinusoidal parallax signature in right ascension to optimize the 
sensitivity of the parallax measurement. The parallax signature 
was followed for two years with a quasi symmetric coverage of 
t=0.0, 0.2, 0.63, 0.64 years with an equal sampling of the min- 
ima and maxima. These two conditions allow one to separate 
the proper motion from the parallax signature. Each observa- 
tion lasted 12 hours and made use of geodetic-like observing 
bloc ks to calibrate the tropospheric zenith delays at each antenna 
(see Reid & Brunthaler"2004t iBrunthaler et"ani2005l: iReid et all 
2009a for a detailed discussion). 

The metha nol masers were first selected from the 

iPestalozzi et al.l ( l2005h database and then observed with the 
Expanded Very Large Array (EVLA, p rogram ABO 16) to ob- 
tain accurate positions (IXuetal.1 12009 b). To find extragalactic 
background sources to serve as reference positions, to verify 
their compact emission, and to obtain obtain positio n with sub- 
arcsecond accuracy, we observed compact NVSS (iCondon et al.l 
1998) sources within 2° of W75N at 5 GHz with the EVN 
in eVLBI mode on December 4, 2008 (program EB039A). 
These observations revealed two compact background sources, 
J2045-H4341 and J2048-H4310 (hereafter J2045 and J2048); we 
also used J 2029-I-4636 (herafter J2029) from the VLBA calibra- 
tor survey (iBeaslev et al.ll2002h separated by 4.°3 from W75N. 
A typical EVN observing run started and ended with a ~1 
hour geodetic-like observing block and about ten minutes for 
fringe-finder observations of 3C 454.3 and J2038-I-5119. The 
remaining time was spent on maser/background-source phase- 
referencing observations. The 6.7 GHz masers in Cygnus X 
and the three background sources were phase-referenced to the 
strongest maser in W 75N, using a switching cycle of 1 .5 min- 
utes. Table[T]lists the source positions, separations from W75N, 
brightnesses, and restoring beam information. DR21 A, DR21B, 
and DR21(OH) are three masers thought to belong to the same 
SFR, but they were observed separately because their separation 
was at the border of the field of view limited by time-smearing 
(~37 arcseconds). 

The observations were performed in dual circular polariza- 
tion and two-bit Nyquist sampling, for an aggregate recording 
rate of 512 Mbps. The data were correlated in two passes at the 
Joint Institute for VLBI in Europe (JIVE) using a one-second av- 
eraging time. The maser data were correlated using one 8 MHz 



' The European VLBI Network is a joint facility of European, 
Chinese, South African, and other radio astronomy institutes funded 
by their national research councils. 
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Table 1. Source parameters of the EVN and the VLBA (AFGL2591) observations 



Source 


R.A. (J2000) 


Dec. (J2000) 


/ 


b 


(f) " 


RA.° 


Brightness* 


Restoring B 


eam* 




(h:m:s) 


(° : ' : ") 


(°) 


(°) 


D 


(°) 


(Jy beam"') 


(mas, mas, deg) 


W75N 


20:38:36.426 


42:37:34.80 


81.87 


0.78 





- 


0.2-23.8 


4.87,3.65,- 


-50.7 


DR21A 


20:39:01.993 


42:24:59.29 


81.75 


0.59 


0.2 


159.5 


0.1 - 1.4 


4.95,3.80,- 


-50.0 


DR21B 


20:39:00.376 


42:24:37.11 


81.74 


0.59 


0.2 


161.3 


0.2 - 1.3 


4.90, 3.77, - 


-52.1 


DR21(OH) 


20:39:01.057 


42:22:49.18 


81.72 


0.57 


0.2 


162.9 


0.02-0.15 


4.90, 3.80, - 


-51.3 


DR 20 


20:37:00.962 


41:34:55.70 


80.86 


0.38 


1.1 


-164.3 


0.04 - 0.8 


4.91, 3.82,- 


-50.8 


IRAS 20290+4052 


20:30:50.673 


41:02:27.54 


79.74 


0.99 


2.1 


-138.0 


0.1 -2.3 


4.91, 3.74, -48.6 


J2045+4341 


20:45:07.055 


43:41:44.51 


83.44 


0.50 


1.6 


48.0 


0.0027 


6.02, 4.47, - 


-39.8 


J2048+4310 


20:48:19.521 


43:10:42.09 


83.41 


-0.28 


1.8 


72.8 


0.011 


4,87, 3,48, -48.0 


J2029+4636 


20:29:18.937 


46:36:02.25 


84.08 


4.48 


4.3 


23.3 


0.041 


4.71, 3.27,- 


-51.4 


AFGL2591 


20:29:24.823 


40:11:19.59 


78.89 


0.71 







17.7 


0.80, 0.40, - 


-14.9 


J2007+4029 


20:07:44.945 


40:29:48.60 


76.82 


4.30 


4.1 


-63.0 


0.992 


0.79, 0.40, - 


-18.6 


J2033+4000 


20:33:03.669 


40:00:24.33 


79.15 


0.04 


0.7 


125.0 


0.025 


0.81, 0.40, - 


-17.9 


J2032+4057 


20:32:25.771 


40:57:27.83 


79.85 


0.70 


1.0 


36.8 








J2033+4040 


20:32:45.342 


40:39:38.12 


79.85 


0.70 


0.8 


53.5 









Notes. Separation, (p, and position angles (east of north), RA. , between the phase reference sources, maser W75 N or AFGL2591, and the 
target. The brightness and restoring beam size (east of north) are listed for the first epoch. 



band with 1024 spectral channels, resulting in a channel separa- 
tion of 7.81 kHz and a channel width of 0.35 km s The back- 
ground sources were correlated in continuum mode with eight 
subbands of 8 MHz bandwidth and a channel width of 0.25 MHz. 
The data were reduced using the NRAO's Astro nomical Image 
Proce ssing System (AIPS) and ParselTongue (iKettenis et al.l 
The data were phase-referenced to the W 75N maser at 
V = 7.1 km s ' and the solutions were transferred to the other 
masers and to the continuum background sources. More de- 
tails o n the reduction and calibration can be found in lRvgl et al.l 
(l20Toh . 

2.2. VLBA water maser observations 

We performed water maser observations using the National 
Radio Astronomy Observatory'^ VLBA in four epochs between 
November 2008 and November 2009 under program BM272H. 
The observation dates were chosen to match the minima and 
maxima of the sinusoidal parallax signature in right ascension 
for optimal paral lax and proper motion analysis, as discussed in 
ISato et al.l ( I2OIOI) . For the VLBA observations, the parallax sam- 
phng was symmetrically covered at t=0, 0.49,0.51, 1.0 years for 
one year. 

To find suitably compact background sources near the maser 
target, we observed a sample of unresolved NVSS sources with 
the VLA in BnA configuration on October 5, 2007, at X and 
K band (program BM272). With these observations, we mea- 
sured source compactness and spectral index and selected the 
strongest, most compact ones to serve as background sources 
(position references) for the parallax observations. Additionally, 
the VLA data were used to obtain subarcsecond accurate posi- 
tions for the background sources and the masers (which were 
also included in these VLA observations). The positions, sepa- 
rations from the maser target, brightnesses, and restoring beam 
information are given in Table [1] 

The VLBA observations, performed at the water maser fre- 
quency of 22.235 GHz, included four geodetic-like blocks for 



^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 



calibrating the zenith delays at each antenna. Two fringe-finders 
(3C 454.3 and 3C345) were observed at the beginning and in 
the middle of each run. The water maser was observed together 
with four background sources, one ICRF calibrator J2007+4029 
dMa et alj ,1998, hereafter J2007), and three quasars selected 
from the VLA observations of NVSS sources: J2033+4000 
(hereafter J2033), J2032+4057, and J2033+4040. We performed 
phase-referencing observations by fast-switching (every 30 sec- 
onds) between the maser and each of the four background 
sources. We used four adjacent subbands of 8 MHz bandwidth in 
dual circular polarization. Each subband was correlated to have 
256 spectral channels, giving a channel width of 0.42 km s"'. 
The data were correlated with the VLBA correlation facility in 
Socorro using an averaging time of 0.9 seconds. The calibration 
was carried out in AIPS following the procedure described in 
Rei d etalJ(l2009ah . 

3. Parallax fitting 

Generally, we detected 6.7 GHz methanol and 22 GHz water 
maser spots in several velocity channels each and toward multi- 
ple locations on the m ap (see Fig.[T]for the methanol masers, and 
see ISanna et aI.ll20Tll for the water masers). All positions were 
determined by fitting a 2-D Gaussian brightness distribution to 
selected regions of the maps using the AIPS task "JMFIT". The 
parallaxes and proper motions were determined from the change 
in the positions of the maser spots relative to the background 
sources (position reference). We fitted the data with a sinusoidal 
parallax signature and a linear proper motion. When the minima 
and maxima of the parallax signature are sampled equally, the 
proper motion and parallax are uncorrelated. Maser spots with 
strongly nonlinear proper motions or a large scatter of position 
about a linear fit were discarded, since these usually reflect spa- 
tial and spectral blending of variable maser features and cannot 
be used for parallax measurements. Only compact maser spots 
with well-behaved residuals were used for the parallax fitting. 

The formal position errors can underestimate the true uncer- 
tainty on the position, since they are only based on the signal- 
to-noise ratios determined from the images and do not include 
possible systematic errors (e.g., from residual delay errors). 
Therefore, to allow for such systematic uncertainties, we added 
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(0,0) = 20:38:36.426, +42:37:34.80 (.12000) _^ (0,0) = 20:39:01.993, +42:24:59.29 (.12000) _^ 
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Fig. 1. 6.7 GHz maser spot maps for W75N, DR21A, DR21B, DR21(OH), DR20, and IRAS 20290 from the first epoch of the 
EVN observations. The spot size is (logarithmicaly) scaled to the peak intensity of the maser spot. 
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error floors in quadrature to the position errors; i.e., we increased 
the positional eiTor for all the epochs by a fixed amount, until re- 
duced values were close to unity for each coordinate. The 
parallax fi tting was then per formed following the same proce- 
dure as in et al.l (l2010i) : 1) we performed single parallax 
fits per maser spot to a background source; 2) we fitted all the 
maser spots (of one maser source) together in a combined par- 
allax fit; 3) when a maser had three or mor e maser spots, we 
perfo r med a fit on the "ave raged data" (see iBartkiewicz et alj 
l2008t iHachisuka et al.l 120091) . However, different maser spots 
can be correlated because an unmodeled atmospheric delay will 
affect all the maser spots of one maser source in the same way. 
Therefore, we multiplied the uncertainty of the combined fit by 
V^, where is the number of maser spots to allow for the like- 
lihood of highly correlated differential positions (procedure 2). 
To obtain the averaged positions for a maser spot (procedure 3), 
we performed parallax fits on all the individual spots and re- 
moved their position offsets and proper motions, after which we 
averaged the positions of each epoch. The last approach has the 
advantage of reducing the random errors, introduced by small 
varia tions in the intern al spot distribution for a maser feature 
(e.g.. lSanna et al.ir2010l) . while leaving the systematic errors un- 
affected. 

The EVN observations used three background sources, of 
which only one, J2045, was usable as a position reference: this 
quasar was close to the maser at the phase reference position 
and unresolved (Fig. |2]i. Quasar J2029 had a too large separa- 
tion from W75N for astrometric purposes (4.°3), which can be 
seen from the near-milli-arcsecond scatter in the variation of the 
relative position versus time of J2009 with respect to J2045 (see 
Fig-EJ. The relative position of J2048 with respect to J2045 (also 
Fig.|3]l also shows a large scatter, thanks to the resolved emission 
of J2048 (Fig.|2]l. Previously, we showed in iRygl et al.l (1201 Ol) 
that at 6.7 GHz some quasars can mimic an apparent motion due 
to changes in the quasar structure. Two or more quasars allow 
one to quantify this effect so that the uncertainty in the maser 
proper motion, which is the error from the derived maser proper 
motion fit and the apparent motion of the quasar, can be esti- 
mated. From the plots of relative positions between J2045 and 
J2029 (J2048 was not considered because it was resolved), we 
estimated the apparent relative motion to be near zero; -0. 1 ± 0.2 
in right ascension and 0.2 ± 0.2 mas yr ' in declination. 

The VLBA 22 GHz observations included four background 
sources, of which two were used for astrometry, J2007 and J2033 
(Fig.|4]i. Two background sources were discarded, one because 
of a non-detection (J203 3 -1-4040) and the other one because of a 
strong structure change during the observations (J2032-1-4057). 
A third source (J2033) was found to be heavily resolved by the 
VLBA (see Fig.|4]l. For this quasar, the position was derived by 
fitting only the peak of the emission in the image, instead of its 
whole spatial structure. 

The parallax and proper motion fitting of the water maser 
was done in the same fashion as for the methanol maser data. 
As the VLBA observations used two background sources, the 
"averaged data" were calculated separately for each background 
source. The combined fit was carried out on the two averaged 
data sets, and the parallax uncertainty was multiplied by Va^ 
for the number of maser spots to account for correlated differ- 
ential positions (see above). In the combined fit, we treated the 
proper motion of the maser separately from the two background 
sources, taking their relative, apparent (linear) motion into ac- 
count (Fig.|5]). This is necessary, because there was a significant 
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Fig. 2. Phase-referenced images of the background sources 
J2045 (top), J2048 (middle), and J2029 (bottom) from the first 
epoch of the EVN observations. The images are in milliarcsec- 
onds offset to the positions in Table[T] The contours start at a 3cr 
level, namely 2.4 x 10""*, 6.6 x 10""*, and 12 x 10""* Jy beam^', re- 
spectively, and increase by V2. The first negative contour (-3cr) 
is shown by dashed contours. The synthesized beam is shown in 
the bottom left corner 5 
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J2045+4341 - J2029+4636 
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Fig. 3. Variation in time of the relative positions of background 
sources J2029 {top panel, red) and J2048 (bottom panel, blue) 
with respect to J2045. The dots mark the right ascension data 
points, while the filled triangles represent the declination data 
points. The solid line shows the right ascension fit, while the 
dashed line represents the declination fit. 



apparent linear motion in the right ascension coordinate (0.32 ± 
0.01 mas yr"'), but none in declination (-0.01 ±0.1 mas yr"'). 

With the parallax and averaged proper motion results we 
calculated the 3D s pace ve locities of the SFRs with respect to 
the Galactic center jReid et al. 2009b). Each SFR has three ve- 
locity vectors: U, the velocity in the direction of the Galactic 
center; V, the velocity in the direction of the Galactic rotation; 
and W, the velocity in the direction of the North Galactic Pole 
(NGP). The maser proper motions that we measure are the ve- 






Fig. 4. Phase-referenced images of the background sources 
J2007 {top) and J2033 {bottom) from the second epoch of the 
VLBA observations. The images are in milliarcseconds offset to 
the positions in Table[T] The contours start at a 3cr level, namely 
5.1 X 10"^ and 1.2 x lO"-' Jybeam"', respectively, and increase 
by V2. The first negative contour (-3cr) is shown by dashed con- 
tours. The synthesized beam is shown in the bottom left corner. 



locity differences between the maser source and the Sun. For 
obtaining the 3D space velocity, we first need to subtract the 
solar motion and then to transform the velocity vector from a 
nonrotating frame to a frame rotating with Galactic rotation. 
In this calculation we used the solar peculiar motion obtained 
by Schonrichet al. (2010) and assumed a flat Galactic rotation 
curve with 6 - 239 km s"' and a solar distance to the Galactic 
center R - 8.3 kpc. These are the rev i sed G alactic rotation pa- 
rameters derived by iBrunthaler et al.l (1201 Ih taking the recent 
revision of the solar peculiar motion into account, from the 
Hipparcos values {Uq, Vq, Wq) = (10.00, 5.25, 7.17) km s"' by 
Dehnen & Binney ( 199^ t o {U^ Vp, Wo)^ (11.10, 12.24,7.25) 
km s ' bv .Schonrich et d] (l2010l) . 
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Table 2. Parallax, proper motions, and space velocities toward Cygnus X 



Source 


Parallax 


Dsun 




<y"<5 > 


I'LSR 


U 


V 


W 




(mas) 


(kpc) 


(mas yr"') 


(mas yr"') 


(km s-') 


(km s"') 


(km s^') 


(kms-') 



W75N 


0.772 


± 0.042 


1 30+""' 


-1.97 


±0.10 


-4.16 ± 


0.15 


9.0' 


-0.5 ± 


1.0 


3.6 


±2.5 


1.2 ± 


0.8 


DR21 


0.666 


± 0.035 


1 50+""** 

-0.07 


-2.84 


±0.15 


-3.80 ± 


0.22 


-3.0' 


-0.0 ± 


1.5 


-8.3 


±2.5 


6.6 ± 


1.3 


DR20 


0.687 


±0.038 


1 46+""'' 


-3.29 


±0.13 


-4.83 ± 


0.26 


-3.0^ 


7.5 ± 


1.5 


-8.8 


±2.5 


5.1 ± 


1.3 


IRAS 20290+4052 


0.737 


± 0.062 


1 36+"'2 


-2.84 


±0.09 


-4.14 ± 


0.54 


-1.43 


1.9 ± 


1.8 


-8.2 


±2.5 


6.0 ± 


2.1 


AFGL2591 


0.300 


±0.010 


Q QQ+0.11 

J.JO g jj 


-1.20 


±0.32 


-4.80 ± 


0.12 


-5.7" 


-12.8 ± 


3.4 


-10.4 


±2.5 


-22.1 ± 


4.3 



References. (D iDickel et an ( ll978h : ISchneider et al.l ( l201(]|) . (2) lMotte et aT] l l2007h . (3) lBronfman et al.Ul996h . (4) lvan der Tak et &\U\99'k . 
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Fig. 5. Variation in time of the relative positions of background 
source J2033 with respect to J2007. The dots mark the right as- 
cension data points, while the filled triangles represent the dec- 
lination data points. The solid line shows the right ascension fit, 
while the dashed line represents the declination fit. 



4. Results 

Table |2] shows the parallax and averaged proper motion results 
for all SFRs and their calculated space velocities, while detailed 
results of the parallax and proper motion fitting are found in 
Table E 



4.1. Cygnus X North: W75N, DR21, andDR20 

W75N has 14 different methanol maser features emitting in an 
local standard of rest velocity (Vlsr) range of 3-10 km s"' with 
32 maser spots in total (see Fig.[T]l. After removing the maser 
spots with nonlinear motions, we were left with ten spots be- 
longing to six maser features for the parallax and proper motion 
fitting. Figure |6] shows the parallax fit of W 75N, resulting in 
0.772 ± 0.042 mas or 1.30+° °^ kpc. While we usually plot the 
parallax fits after removing the proper motion, in this figure we 
also show a parallax fit, to one of the maser spots of W 75N, 
without the removal of the proper motion. Since different maser 
spots can have different proper motions, it is not instructive to 
plot the parallax fits to all the maser spots when including the 
proper motion. 

For DR21A we found 18 maser spots belonging to seven 
maser features, covering an LSR velocity range of -9.4 to 



J2045+4341 



J2045+4341 



O 





2009 2010 

Time (years) 



2009 



2010 
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Fig. 8. Results of the parallax fit for DR 20 after removing proper 
motions and positional offsets. The declination data has been off- 
set for clarity. Left panel: Combined fit on 3 maser spots with 
respect to J2045. Right panel: Fit to the averaged data of the 3 
maser spots. The dots mark the data points in right ascension, 
while the filled triangles mark the declination data points. The 
solid lines show the resulting parallax fit in right ascension and 
the dashed lines show the fit in declination. 



-5.2 km s"' (see Fig.[T]). The parallax fitting was based on 
seven spots in four maser features. Toward DR21B we found 
seven spots in three maser features with Vlsr between 3.3 and 
5.0 km s ' (see Fig.[TJ, but only one maser spot had a good paral- 
lax fit. We also detected methanol masers in DR21(OH) where 
we used three maser spots with Vlsr=-2.5 to -3.5 km s"' for 
fitting a parallax signature. If we assume that the masers in 
DR21A, DR21B, and DR21(OH) belong to the same SFR - 
which is likely since their parallaxes are consistent: 0.686 ± 
0.060 for DR21A, 0.705 ± 0.072 for DR21B, and 0.622 ± 
0.055 for DR21(OH) - the resulting parallax fit for DR21 
becomes 0.666 + 0.035 mas, or 1.50^1}"^ kpc. The connection 
between DR2I and P R 21 (OH) has already been noted by 
[Schneider et al ] (120061) . since these SFRs are part of the same 
elongated filament seen in CO emission. The parallax fits to the 
combined and averaged data for DR 21 A and DR 21(OH), the fit 
to the individual spot of DR21B, and their overall combination 
(i.e., for the DR 21 SFR) are shown in Fig.|7] 

Six maser spots in the velocity range of -4.8 to -2.0 km s"', 
belonging to two maser features that were found in DR 20 (see 
Fig.lTJ. The parallax fit was based on one maser feature with 
three maser spots. We find a parallax of 0.687 ± 0.038 mas, cor- 
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Fig. 6. Results of the parallax fit for W 75N. Left panel: Results of the parallax fit for the maser spot at Vlsr - 9.56 km s ' without 
the removal of the proper motions. The right ascension and declination data has been offset for clarity Middle and right panel: 
Results of the parallax fit for W75N after removing proper motions and positional offsets. The declination data has been offset for 
clarity. Combined fit on 10 maser spots with respect to J2045 {middle panel). Fit to the averaged data of the 10 maser spots {right 
panel). The dots mark the data points in right ascension, while the filled triangles mark the declination data points. The solid lines 
show the resulting parallax fit in right ascension and the dashed lines show the fit in declination. The scale of the Y-axis of the 
middle and right panels is different from the scale of the left panel. 



responding to a distance of 1 .46^° kpc (see the parallax fit in 
Fig-EJ- The second epoch was omitted from the parallax fitting, 
because the maser spots in DR 20 had an outlying data point in 
proper motion caused by large residuals in atmospheric delay in 
that epoch's data. 



4.2. Cygnus X South: AFGL2591 and IRAS 20290+4052 

The parallax of AFGL2591 was based on VLBA data us- 
ing 22 GHz water masers. AFGL2591 has a very rich water- 
maser spectrum of 80 maser featur es with Vlsr from -34 to 
-0.4 km s ' (see San na et alj [201 iV The parallax and proper 
motion fits were based on six maser features and resulted in a 
parallax of 0.300 ± O.OlOmas or 3.33+° ][ kpc. Figure |9] shows 
the parallax and proper motion fit of AFGL2591 with respect 
to the two background sources. The parallax fit to each back- 
ground source separately was 0.302 ± 0.009 mas (J2007) and 
0.299 + 0.002 mas (J2033). The difference between the two par- 
allax fits is due to the apparent motion observed between the 
two background sources. The effect on the parallax can be ap- 
proximated by the difference between the two results, namely 
0.003 mas. The previous distance measurements of AFGL2591 
were near 1.6 kpc. Assuming tha t the SFR is associated to 
IC 1318c (IWendker & Baarslll974l) . of which the distance was 
determined bv lDickel et al.l (1969). AFGL2591 was thought to 
be at 1.5 kpc. Alternatively, Dame & Thaddeus (1985) suggest 
a distance of 1.7 kpc based kinematic distances to CO clouds. 
However, one can find a m uch wider spre ad of distances from 
1 to 2 kpc in the literature (Ivan der Tak et al. 1999). The water 
maser parallax puts AFGL2591 a factor 2.22 (assuming 1.5 kpc) 
farther away, at 3.33 kpc, which implies a dramati c change for 
the p hysical properties of the SFR (discussed in I Sanna et"al] 
120111) . 

Finally, for IRAS 20290, we found two methanol maser fea- 
tures composed of a total of eight spots with Vlsr from -6.2 to 
-3.8 km s ' (see Fig.[T]i. Due to large residuals for the fit, only 
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Fig. 9. Results of the 22 GHz water maser parallax fit for 
AFGL 259 1 (VLBA data) after removing the proper motions and 
positional offsets. The declination data has been offset for clarity. 
Left panel: Combined fit on 6 maser spots with respect to J2007 
and J2033. Right panel: Fit to the averaged data of the 6 maser 
spots. The dots mark the data points in right ascension, while the 
filled triangles mark the declination data points. The solid lines 
show the resulting parallax fit in right ascension and the dashed 
lines show the fit in declination. 



two maser spots, belonging to the same feature, were suitable 
for parallax fitting resulting in a parallax of 0.737 ± 0.062 mas or 
1.36+Q [j kpc. The parallax fit is shown in Fig.fTOl 
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Fig. 7. Results of the parallax fit for DR21 after removing proper motions and positional offsets. The declination data has been 
offset for clarity. 7b;? row, from left to right, first panel: Combined fit on 7 maser spots from DR 21 A with respect to J2045. Second 
panel: Fit to the averaged data of the 7 maser spots of DR21A. Third panel: Single fit to DR21B maser spot at 3.96 km s Bottom 
row, from left to right, first panel: Combined fit on 3 maser spots in DR21(0H). Second panel: Fit to the averaged data of the 3 
maser spots of DR21(OH). Third panel: Combined fit on the averaged data set of DR21A, DR21(OH), and the single maser spot 
of DR21B. The dots mark the data points in right ascension, while the filled triangles mark the declination data points. The solid 
lines show the resulting parallax fit in right ascension and the dashed lines show the fit in declination. 
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Fig. 10. Results of the parallax fit for IRAS 20290 after removing 
proper motions and positional offsets. The declination data has 
been offset for clarity. Shown is the combined fit on 2 maser 
spots with respect to J2045. The dots mark the data points in 
right ascension, while the filled triangles mark the declination 
data points. The solid lines show the resulting parallax fit in right 
ascension and the dashed lines show the fit in declination. 
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Table 3. Detailed results of parallax and proper motions measurements 



Maser Background 


vlsr 


Parallax 






source 


(km s'^) 


(mas) 


(mas yr""') 


(mas yr"') 


W75N J2045+4341 


9.56 


0.750 ± 0.070 


-2.01+0.13 


-4.33 + 0.10 




9.23 


0.800 ± 0.067 


-1.99 ±0.11 


-4.37 + 0.11 




9.23 


0.730 ± 0.079 


-2.03 ±0.13 


-4.37 + 0.15 




9.23 


0.748 ± 0.062 


-1.85 ±0.09 


-4.09 ±0.19 




8.88 


0.879 ± 0.060 


-1.97 ±0.09 


-4.42 ±0.13 




8.88 


0.733 ± 0.068 


-1.87±0.11 


-4.11 ±0.13 




6.77 


0.764 ± 0.068 


1 1 O 1 A 1 A 

—2.18 ± O.iU 


A i n 1 n OA 

-4.1 / ± U.zU 




/111 
4.31 


\jJ25 ± 0.072 


1 AO 1 All 

— 1.99 ± 0.11 


1 AT 1 AIT 

-3.93 ± 0.17 




3.26 


A 'lAdZ 1 A AC/C 

0.746 ± 0.056 


1 OA 1 A nA 

-1.89 ± 0.04 


'J (\A 1 A 1 

—3.94 + 0.12 




2.91) 


A OA1 1 A nCdZ 

O.oOl ± 0.056 


1 A^C 1 A AO 

-1.96 ± 0.08 


o oc 1 n tn 

-3.85 + 0.17 




Combined fit 


U.//Z + U.Uo j 








Averaging data 


0.772 ± U.U42 






< n> 






— 1.9/ ± O.iU 


4 t /Z 1 A 1 C 

—4.10 + 0.15 


DR21A J2045+4341 


-5.54 


0.751 ±0.094 


-2.67 ±0.15 


-3.22 + 0.18 




-5.89 


0.706 ± 0.087 


-2.76 + 0.14 


-3.12 + 0.15 




-6.59 


0.709 ± 0.094 


-2.94 + 0.15 


-2.94 + 0.24 






718-1-0 nfis 

U. / lo ± U.UOo 


9 07 -1- n 1 1 

— Z.O / ± U. 1 1 


9 QQ -t- n 17 




-9.05 


0.710 ± 0.11 


-2.59 + 0.17 


-3.87 + 0.38 




-9.05 


0.719 ± 0.11 


-2.78 ± 0.17 


-3.67 ± 0.35 




-9.40 


0.600 ± 0.088 


-3.26 ± 0.14 


-4.37 ± 0.14 




Combined fit 


713 + 091 








Averaging data 


0.686 + 0.060 






< u > 






-2.84 ±0.15 


-3.56 ± 0.25 


r^T? '>1R Ton/K_i_/n/i 1 

UK. Zlo JZU4j+4o41 


j.yO 


A TrtC _i_ A AT'! 
U. /Xfj ± U.U/Z 


T QQ -.All 
— Z.OO ± U. i 1 


— 4.:>0 ± V.Lj 


DR21(OH) J2045+4341 


-2.72 


0.694 ± 0.100 


-2.31 ± 0.16 


/I A 1 /' 

-4.29 ± 0.16 




-3.07 


ACTI 1 A 111 

0.533 ±0.111 


-3.11 ± 0.19 


-4.51 ± 0.17 




-3.43 


0.572 ± 0.104 


-3.05 ± 0.16 


-4.51 ± 0.22 




Combined fit 


A cat 1 A 1 nA 

0.591 ± 0.104 








Averaging data 


U.022 ± U.055 










n tmti n n'\^ 

U.OOO ± 






< yCi > 






-2.84 ±0.15 


-3.80 ±0.22 


UK2\) Jz(J43+4341 


-3.78 


0.634 ± 0.044 


-3.16 ± 0.06 


A 1 A OA 

-4.82 ± 0.39 




-4.13 


A TIA 1 A ACiC 

0.729 ± 0.056 


At 1 A AO 

-3.41 + 0.08 


A OC 1 A A^ 

—4.85 + 0.42 




A AH 


U. / ZD ± U. 1 jO 


"J 70 -1- n 91 
— j.zy ± U.Zl 


4 S 1 J- n 4^ 




Combined fit 


0.700 ± 0.085 








Averaging data 


0.687 ± 0.038 






< yU > 






-3.29 ± 0.13 


-4.83 + 0.26 


IKAo 202!*0+4052 J 2045 +4341 


C CA 

-5.54 


A /;A 1 1 A AjCA 

0.691 ± 0.060 


O 1 A A A 

-2.82 ± 0.09 


/I 1 T 1 A CA 

—4.1 / + 0.50 




on 

-5.89 


A TOC 1 A AiCO 

0.785 ± 0.063 


O C 1 A AA 

-2.85 ± 0.09 


-4.11 + 0.57 




Combined fit 


A TIT 1 A AiCI 

0.7 j7 ± U.Ud2 






< yCi > 






T o /I 1 n nn 
-2.o4 ± 0.09 


/I 1 /I 1 A C/1 

—4.14 + 0.54 


AFGL2591 J2007+4029 


-4.58 


0.267 ± 0.009 


-0.656 ± 0.025 


-5.906 + 0.051 




-5.00 


0.308 ±0.014 


-1.083 ±0.037 


-4.019 + 0.169 




-5.84 


0.290 ± 0.007 


-1.271 ±0.018 


-4.152 + 0.169 




-7.53 


0.317 ± 0.004 


-0.955 ±0.011 


-4.715 + 0.071 




-7.95 


0.280 ± 0.005 


-1.385 + 0.015 


-5.745 + 0.155 




-10.90 


0.348 ± 0.008 


-0.884 + 0.021 


-4.398 + 0.113 




Averaging data 


0.302 ± 0.009 


- 


- 


J2033+4000 


-4.58 


0.263 ± 0.003 


-0.980 ± 0.009 


-5.863 + 0.169 




-5.00 


0.301 ±0.009 


-1.407 ±0.025 


-3.976 + 0.025 




-5.84 


0.287 ±0.001 


-1.595 ±0.002 


-4.109 + 0.028 




-7.53 


0.315 ±0.001 


-1.279 ±0.004 


-4.672 + 0.141 




-7.95 


0.277 ± 0.001 


-1.709 ±0.001 


-5.702 + 0.031 




-10.90 


0.345 + 0.002 


-1.208 ±0.005 


-4.355 + 0.078 




Averaging data 


0.299 + 0.002 






Botii QSOs 


Combined fit 


0.300 + 0.010 






<//> 






-1.20 + 0.32 


-4.80 + 0.12 



Notes. The resulting parallax fit for DR21 is based on the averaged data of DR21A, DR21(0H), and the single maser spot of DR21B. 
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5. Is Cygnus X one region? 

AFGL2591, supposedly located in the Cygnus X South region, 
has a much farther distance than previously assumed, which 
implies that AFGL2591 is not a part of a single Cygnus X 
complex. AFGL2591 is perhaps part of the Local Arm, which 
would then extend to greater distances from the Sun than cur- 
rently thought. While we need more data to investigate this issue 
with stronger statistical support, we note that AFGL 2591 is not 
the only source found in this region (see Fig.fTTI). Two recent 
parall ax measurements, namely ON2 at 3.83 kpc dAndo et alj 
1201 Ih and G75.76+0.35 at 3.37 kpc (Xu, private comm.; see 
also the Bar and Spiral Structure Legacy survey, BeSSeL, web- 
sit^, locate these SFRs close in space to AFGL 2591 (within 
3° on the sky). AFGL 2591 has a Vlsr similar to other SFRs in 
Cygnus X North, even though it is projected against the Cygnus 
complex and not part of it. The space motions of AFGL 2591, 
though, are very different from the sources in Cygnus X as can 
be seen in Fig. [12] We note that the distance measurement of 
AFGL 2591 is not affected by using different VLBl arrays and 
maser transitions than for the other SFRs. Maser parallaxes, both 
the water and the 6.7 (and 12.2) GHz methanol masers have 
been shown to produce robust distance measurements: for ex- 
ample, the distance to W3 (OH), measured wit h both 12.2 GHz 
methanol (IXu et al.ll2006l) and 22 GHz water (iHachisuka et al.1 
I2OO6I) masers, yi elded the same n umber, whereas the 6.7 GHz 
maser distances (Rygl et al."2010) also agreed with the VERA 
water maser distances (Sato et al. 2008; Nagavama et al. 201 1). 

The most important result of this study is that Cygnus X 
North is one physically related complex of SFRs, including 
W75N, DR21, DR20, and IRAS 20290 (and therefore proba- 
bly also the Cygnus OB 2 association), located at 1.40^QQgkpc. 
This is an average of the individual distances, which range from 
1.30 to 1.50 kpc. Our data are consistent with a single distance 
for these sources, within measurement uncertainty. We note that 
our distance to the Cygnus X complex is similar to the photo- 
metric distance of 1.5 kpc obtained b v iHansonl (12003 ) toward the 
Cyg OB 2 association. Compared to the extent on the sky of the 
SFRs mentioned above, 25 x 60 pc, the distance spread of 200 pc 
is a factor ~3.5 wider (however, the measurement uncertainty is 
a factor 1-2 times the angular extent). The parallax results can 
possibly be extend ed based on cont i ngent spatial-velocity struc- 
tures identified by [Schneider et al.l (l2006l) in the CO emission. 
Following these authors DR22, DR23, DR 17, and AFGL 2620 
(their CO groups 1 and II) should also be part of Cygnus X North 
and thus at the same distance. 

Our results did not find any evidence of a southern coun- 
terpart to Cygnus X North, since AFGL 2591 was found to be 
much more distant and other SFRs toward Cygnus X South 
(except for IRAS 20290 at Z=79.°7, which was found to be at 
the same distance as Cygnus X North) were not included in 
this work. [Schneider et al.l ([2006h find that most of the mass in 
Cygnus X South is contained in a group of molecular clouds: 
DR4, DR5, DR12, DR13, and DR 15 (their CO group IV). 
Parallax measurements to one of these SFRs then could con- 
firm that Cygnus X South is connected to Cygnus X North, thus 
forming one of the largest known giant molecular cloud struc- 
tures in the Milky Way of ~ 2.7 x 10^ Mq dSchneider et al.|[2006[ 
rescaled to 1.4 kpc). 

The distances of the SFRs in Cygnus X fit well with the 
trajectory of the Local Arm between 0.5 and 2.5 kpc in the 
X coordinate (see Fig.[TT[|. defined by measurements to CepA 
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Fig. 11. An artist's impression of the Galactic plane (im- 
age credit: R. Hurt, NASA/JPL-Caltech/SSC) as seen from 
the North Galactic Pole with the Galactic center (not shown) 
located at (0,0) kpc. The Sun is located at (0, 8.3) kpc 
and marked by a black circle around a black dot. Dots 
mark the sources to which an accurate distance is known 
in the Local Arm: green dots mark sources from this 
work; red dots mark sourc e s from the Uteratu r e (from left 
to right): 'Choi et al. (2008), Reidetal. (2009a), Hirota et a^ 
((2007); Menten et al. (2007), .Hi rota et al. ( 2008a) . ,Hirota eta" 
(201 1), Moscadel U eit alj (|2009iriR vgl et alj| 
(2008b), Reid et al.'('201ll). [Xuet al.i ' (i2009ah . [Rvgl et al 
Nagavama et al. (2011), Xu, private comm., and ,Ando et akl 
(1201 II) . For a comparison we also plot the sources in 



Hirota et al 
.' (201(1: 



the Perseus Arm (fr o m left to right) : [Reid et alJ ([2009a|) . 



iNiinumaetd] ([201 Ih . [Oh et al 
Hachisuka et al.'('2006'),'Asaki et al 
(2007), Sato et al. (2008); Rygl et al 



j2om 

(2010), 



IXu et al.l (I2OO6 ) 



Moellenbrock et al 
(2010), Mos cadelH et 4 
( I2009D . and lOhetal] ( [2010) . the Local and Perseus Arms are 
indicated. 



(iMoscadeUi et alJ E009l) . L1206 (iRvgl et all l20Toh . L 1448C 
( [Hirotaetal.[[2011 [) , Cygnus X-1 (lo cated ~ 7° away from the 
Cv gnus X region. [Reid et all [201 ll) . ONI dRvgl et alJ [20101: 
Nag avama etalJl20 111) . andG59.7-H01 (IXu et al.ll2009al) . 

There are two explanations proposed for a single con- 
nected Cygnus X region. First, a superbubble driven by the fa- 
mous Cygnus Loop sup ernova remnant (e.g., [Walsh & Brownl 



1955 : Cashetal.1 19801) . which was recently discarded by 



Uvaniker et alJ ([20011). and, second, an expanding Stromgren 



http://www.mpifr-bonn.mpg.de/stafF/abrunthaler/BeSSeL/index.shtml 



sphere (McCutcheon & Shuteii [1970). Our Cygnus X sources 
(all sources except AFGL 2591) have overlapping distance un- 
certainties, so we cannot use their distances for studying the 
structure of the complex. However, the UVW space motions of 
the sources projected on the Galactic plane give an impression 
of the dynamics of the Cygnus X complex. We plot the result- 
ing UVW space motions toward Cygnus X in Fig. [T2| Apart 
from a clearly different behavior of AFGL 2591, we find that 
W75N and DR21 have a dominant motion toward the NGP, 
while DR20 and IRAS 20290 are moving toward the NGP and 
the Galactic center These proper motions of the Cygnus X SFRs 
do not suggest a common expansion center, which an expanding 
Stromgren sphere should have, so more data are needed to un- 
derstand the formation of the Cygnus X region. 

Additionally, we found that the space velocity V, in the di- 
rection of the Galactic rotation (after subtracting the Galactic ro- 
tation of y = 239km s"'), lies between -10.5 and -8.0km s"' 
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Fig. 12. Midcourse Space experiment (MSX) 8 /zm image of tlie Cygnus X region overlaid with the resulting UVW space motions 
for each source. The white triangles mark the water maser (AFGL2591) and the methanol maser sources, with their distances 
indicated in kpc. 



for all sources except W75N, where V - +3.6 km s" . It seems 
that most of the star-forming gas is moving with the same 
Galactic orbital velocity, lagging some 9 km s"' behind circular 
orbits, as found b y most parallax studies of massive SFRs (e.g., 
iReid et al."2009b') after taking the revised Solar motion into ac- 
count (Schonrich et al. 2010,) . 

6. Summary 

We measured the trigonometric parallaxes and proper motions 
of five massive SFRs toward the Cygnus X star-forming com- 
plex using 6.7 GHz methanol and a 22 GHz water maser. 
We report the following distances: 1.30;^qq7 kpc for W75N, 
1.46!|J{J^kpc for DR20, 1 .50+JJ JJ^ kpc for DR21, 1.36+° [2kpc 
for IRAS 20290+4052, and 3.33+° [] kpc for AFGL 2591. While 
the distances of W75N, DR20, DR21, and IRAS 20290+4052 
are consistent with a single distance of 1.40 + 0.08 kpc for the 
Cygnus X complex, AFGL 2591 is located at a much greater dis- 
tance than previously assumed. The space velocities of the SFRs 
in Cygnus X do not suggest an expanding Stromgren sphere. 
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